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Abstract 

We consider an optimal control problem for the obstacle problem with an elliptic variational 
inequality. The obstacle function which is the control function is assumed in H 2 . We use an ap- 
proximate technique to introduce a family of problems governed by variational equations. We prove 
optimal solutions existence and give necessary optimality conditions. 
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1 Introduction 

The study of variational inequalities and free boundary problems finds application in a variety of dis- 
ciplines including physics, engineering, and economies as well as potential theory and geometry. In the 
past years, the optimal control of variational inequalities has been studied by many authors with differ- 
ent formulations. For example optimal control problems for obstacle problems (where the obstacle is a 
given (fixed) function) were considered with the control variables in the variational inequality. Roughly 
speaking, the control is different from the obstacle, see for example works by [4], [8], [15], [16] and the 
references therein. 

Here we deal with the obstacle as the control function. This kind of problem appears in shape optimization 
for example. It may concern a dam optimal shape. The obstacle gives the form to be designed such that 
the pressure of the fluid inside the dam is close to a desired value. This is equivalent in some sense to 
controlling the free boundary [10] . 

The main difficulty of this type of problem comes from the fact that the mapping between the control 
and the state (control-to-state operator) is not differentiable but only Lipschitz-continuous and so it is 
not easy to get first order optimality conditions. 

These problems have been considered from the theoretical and/or numerical points of view by many 
authors (see for example Adams and Lenhart [3], Ito and Kunisch |12j). They have used either an ap- 
proximation of the variational inequality by penalization-regularization or a complementarity constraint 
formulation. Adams et Lenhart [3] consider optimal control problem governed by a linear elliptic varia- 
tional inequality without source terms. The main result is that any optimal pair must satisfy "state = 
obstacle". Adams and Lenhart [2] treat control of H 1 — obstacle, and convergence results in [2], [3] are 
given under implicit monotonicity assumptions. 

Ito and Kunisch [12] consider the optimal control problem to minimize a functional involving the H 1 
norm of the obstacle, subject to a variational inequality of the type y 6 argmin{a(z) — (/, z)\z < ip} 
in a Hilbert lattice H. Under appropriate conditions, they show that the variational inequality can be 
expressed by the system Ay + A = /, A := max(0, A + c(y — -0))- Smoothing the max-operation, this 
system is approximated by a semilinear elliptic equation containing only smooth expressions. Passing to 
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the limit, the optimality system of the associated differentiable optimal control problem is used to derive 
an optimality system of the original nonsmooth control problem with only H ^regularity for the obstacle. 
Bergounioux and Lenhart [6j, [7] have studied obstacle optimal control for semilinear and bilateral obstacle 
problem, where the admissible controls (obstacles) are i/ 2 -bounded and the convergence results are given 
with a compactness assumption. Yuquan and Chen [17], consider an obstacle control problem in a 
elliptic variational inequality without source terms. We can see also quote the paper of Lou [H] for more 
generalized regularity results. In this paper we consider an optimal control problem: we seek an optimal 
pair of optimal solution (state, control), when the state is close to a desired target profile and satisfies 
an unilateral variational inequality with a source term, and the control function is the lower obstacle. 
Convergence results are proved with compactness techniques. 

The new feature in this paper is the regularity on the control function (obstacle) and an optimality 
conditions system more complete than the one given in |17| . 

Let us give the outline of the paper. Next section, is devoted to the formulation of the optimal control 
problem, we give assumptions for the state equation, and give preliminaries results. In section 3 we 
study the variational inequality, give control-to-state operator properties and assert an existence result 
for optimal solution. The last section is devoted to optimality condition system. 



2 Optimal control problem 

Let SI be an open bounded set in M" (n < 3), with lipschitz boundary dfl. We adopt the standard 
notation H m (Q) for the Sobolev space of order minO with norm ||-||jym(n), where 

H m (Sl) := {v | v e L 2 (Sl), d q v e L 2 (Sl) Vq, \q\ < to} , 

and 



H?(il) := \v | v G H™(Sl), ;1 , 



drj k 



= 0, <k < to - 1 

an 



defined as the closure of ^ (SI) in the space H m (SI), where @(SY), the space of C°° (Sl)-functions, with 

compact support in SI (see for example [1]). We shall denote by \\-\\ v , the Banach space V norm, and 

l/p 



ljy(n) 



the p — summable functions u : £1 — > R endowed with the norm ||u||ip(n) := (^J \u (x) \ p dx 



for 1 < p < oo and ||w||l°o(o) : — esssup|M(s) | for p — oo. In the same way, (•,•} denotes the duality 

product between H~ 1 (Sl) and Hq(SY), and (•, •) the L P (ST) inner product. It is well known that Hq(Q) <^-* 
L 2 (Sl) <—> H^ 1 (Sl) with compact and dense injection. We consider the bilinear form a(-,-) defined on 
H*(Sl) x ffi(fi) by 

a(u,v):= J2 / a ^JfT-^r dx + J2 a ^ v dx + / a uv dx, (2.1) 



n dxi dxj Jq dxi j n 



where 



a ,ai,aij G L°° (SI) 



V <,,,(),(), ■ m V Of. to>0, a.e. infi, M6 e 



(H) 



Moreover, we suppose that aij G C 01 (S),) (the space of lipschitz continuous functions in SI, where fl is 
the closure of SI) and that do is nonnegative to ensure a good regularity of the solution (see for example 
[13J). We suppose that the bilinear form a(-, ■) is continuous on Hq(Q) x Hq(SI) 

3M > 0,V^,^ G Ht(n),\a(<p,iP)\ < M M Hl{n) M H i {n) , (2.2) 
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and coercive on i?o(Q) X #o(Sl) 

3m > 0,^ e H^n),a^,^) > mM 2 Hlom . (2.3) 

We call A G £(iJ* (SI), i/^^S!)) the linear (elliptic) operator associated to a(-,-) such that (Au,v) := 
a(u,v). We note that the coercivity assumption (|2,3p on a implies that 

v>Gi?o(0), {Atp,ip) >m\\tp\\ 2 H i {n) . 

For any ip G Hq (SI), we define 

K(<p) := {y G H^(Q) \y>p a.e. in Si} , 
and consider the following variational inequality 

a (y, v - y) >(/,«- y) , -iv G Kbp), (2.4) 

where / belongs to L 2 (SI) as a source term. In addition c or C denotes a general positive constant 
independent of 8. 

Theorem 2.1. Under the hypothesis \2.2\ and l|2.3j) . /or any / € L 2 (J!) and 99 € i?Q (SI), i/ie variational 
inequality l|2.4p . /ias a unique solution y in IC (ip). In addition if tp belongs to H 2 (SI), i/ie solution y belongs 
to H 2 (Q)nH^{fl). 

Proof. See [9]. □ 

From now we define the operator T (control-to-state) from H 2 {Q) n (SI) to if 2 (SI) n i/g(f2), such that 
y '■= T (tp) is the unique solution to the variational inequality (|2.4p . 
Now, we consider the optimal control problem (T 3 ), defined as follows 

minj.%) := \ jjp {<p) - zf + % ^(A^) 2 dx^j , p e W arf | , (P) 

where v is a given positive constant, z £ L 2 (SI) and £/ a d (the set of admissible control) is a closed convex 
subset of i/ 2 (SI) n -ffg(Sl): we seek an obstacle (optimal control) cp* in Had, such that the corresponding 
state is close to a target profile z. In the sequel we set U := H 2 (Q) n Hq(S1). 

3 Approximation of problem (V) 

3.1 Approximation of operator T 

The obstacle problem (|2,4|1 can be equivalently written as follows 

Ay + dI K{v) (y) 9 / in SI, y = on 5S1, (3.1) 

where 

dl/c( v ) (y) = dI K + iy _ v) (y) := {v G £ 2 (SI) | v G /3 D (y - <p) , a.e. in Si} , 

and 

/C+ := {y G Ho(Sl) I y > 0, a.e. in Si} , 
and f3 : M. — >2 R is the maximal monotone (multivalued) graph, 
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{0 if r > 
R- if r = 
if r < 0. 



Equation l|3.ip . can be approximated by the following smooth semilinear equation 



Ay + p s (y - cp) = f in O, y = on SO, 
where /3s is an approximation of /3 Q . One possible approximation of (3 si given as follow 



(3.2) 




if r > 

-r 2 if r £ [-|,0] 
r + | ifr<-|, 



Where <5 > and we note that r s := | min{0, r} < (r) < and /3 6 < ^°° (E) and /3£ is given by 



As l3$(- — tp) is nondecreasing, it is well known (see for example [H]), that the boundary value problem 
(pT2|) admits a unique solution y s in H 2 (ty n #o( fi ) for a fixed f m # 2 (^) H #o(fi) and / in L 2 In 
the sequel, we set y 5 := T 5 (y). We recall the following continuity results [7] 

Theorem 3.1. i<br ant/ pair (j/i, ipi) inU xU, that satisfies l|3.2p where i = 1, 2. We </e£ 

||S/2 ~ 2/l|liJi(Q) < i«5 \\<f2 - <Pl\\ L 2 (Q) , 

where Ls ■= maxjl, and m is i/ie coercivity constant of a. 
Proof. From (|3.4ap . we obtain 

a {yi -yi,v) + (Ps (2/2 - <P2) - As (2/1 - Vi) , v) = 0, Vw£ ZY. 
with v = 2/2 — 2/1 , we write 

a (2/2 - 2/1,2/2 - 2/1) + (#5 (2/2 - ^2) - As (2/1 - ^2) ,2/2 - 2/1) = 0. 
Since /3a is nondecreasing, by the hypothesis 112.21) and l|2.3p . we deduce 

\\V2 - 2/i|| ff i ( n) < L s Wfz ~ <Pi\\ L *(n) » 
where La := max {l, □ 

Theorem 3.2. Let ip s in H 2 (Q) n i?o(il) 6e a strongly convergent sequence in Hq (ft) to some <p as 5 
tends to 0. Then the sequence y s :— T s (<p 5 ) strongly converges to y :—T ((p) in Hq (f2). 

Proof. For every ip s in H 2 (fl) n H^(fl), we set y s := T 5 (<//), then for any 2/ in Hq (fi) the equation 
p.2p is equivalent to 





•)=(/,«), VveH^n). 



(3.3) 



a (2/, / - y s ) + [ As (y s - <p s ) (<p s - y s ) dx 



[ f(p 5 -y 5 ) dx 
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We know by the definition of ps, thats if y s (x) — (p 5 (x) > 0, we have 

p s (y s (x) - <p s (x)) = 0, 
otherwise, we get P$ (y 5 (x) — p s (x)) < 0. Then we deduce that in all cases, we have 

p s (y 5 - </) (y s ~ <p 5 ) > a.e. in Si, 

that yields 

a(y s ,y s ) <a(y s ,p s ) + (f,y 5 -p s ), 
with the hypothesis (12.21) and j|2.3|) . we deduce (estimation regularity) 

\\y 5 \\HHn)^ C h S \\m(n)> ( 3 - 4 ) 

where C is a constant only depending on / and a. We know that if <p s is strongly convergent in Hq (Si) 
then ip 5 is bounded in Hq (Si), and by l|3.4p we deduce that y s is convergent to some y as 6 tends to 
weakly in Hq (Si) and strongly in L 2 (Si), 

Let w in K (p), and choose v s = max (w, i/? 5 ) . We have v 5 in JC (<p s ) and that w" 5 is convergent to v strongly 
in Hq (SI). Equation l|3.3p with v — v s — y 5 gives 

a (y s , v s - y 5 ) + [ p & (y s - ip 5 ) (v 5 - y s ) dx = [ f (v 5 - y 5 ) dx. 
Jn Jn 

• If y 5 < <P S , therefore Ps (y s - p 5 ) < and (v s - y 5 ) > 0, we deduces that Ps (y s - <p s ) (v 5 - y 5 ) < 
0. 

• If y S > V s , then y s - p 5 > 0, therefore Ps (y 5 - <p s ) = 0. 

So we deduce that in all cases we have Ps (y s — (p s ) (v 5 — y s ) < 0, and we get 

a(y 5 ,y 5 )<a(y 5 ,v 5 ) - (f,v 5 ~y s ). 
Passing to the limit and using the lower semi-continuity of a gives 

u yy, y ) n luuuu u y[' 5 " ! 1 1 " ' " 

and 



a(y,y) < liminf a(y s ,y s ) < Uminf a (y S ,v s ) - (f,v s - y A ) =a(y,v) -(f,v-y). 

o— >0 ()—*{) 



a(y,v-y)>(f,v-y), W e fC (p) . 

It remains to prove that y 5 tends to y, strongly in Hq (Si). By using the fact that w s — m&x(y,ip s ) 
converge to y strongly in Hq (Si) it is sufficient to prove that w s — y s converge to strongly in Hq (St). 
From equation l|3.3p we get 

f 5 S S 8\ ( 5 5 5\ ( 5 5 5\ 

a(w - y , w -y ) = a(w ,w - y ) - a(y ,w - y ) 

= a (w s ,w s -y s )+ f p s (y s - p s ) (w s - y s ) dx - [ f (w s - y s ) dx. 
Jn Jn 

As previously we deduce that 

a (w s - w 5 ~y 5 )<a (w 5 ,w s - y s ) - (/, w 5 - y s ) , 



from the hypothesis (|2 .3[) . we get 

II 5 5 1| 2 ^ ( 5 5 5\ I r 5 S\ 

m \\ w ~y \\m(n) < a (w ,w -y ) - (f,w -y ) 



□ 
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As a consequence of the previous theorem, we obtain the following corollaries 
Corollary 3.1. For any tp s in H 2 (Si) n H% (Si), y s := T s (p s ) belongs to H 2 (St) n H% (SI). 
Proof. Since fa {v & - P & ) and / belongs to L 2 (SI), then Ay 5 G L 2 (SI) and y 5 G i? 2 (fi). □ 
Corollary 3.2. For any <p inU a d, the sequence y s ■= T s (ip), converges toy :— T (ip), strongly in Hq (St). 
Corollary 3.3. There exists a constant C depending only on f and a, such that for any p in U, we get 

\\^m H ^)<C{aJ)\W\\ Hlm . 

Proof. We choose p 5 = (fi, and y s := T 5 (ip), as we know that y s converges to T (ip) strongly in Hq (St), 
we pass to the limit in p,4p . □ 



Theorem 3.3. T is continuous from U endowed with the sequential weak topology of H 2 (SI) to Hq (SI) 
endowed with the sequential weak topology. 

Proof. Let p k be a sequence that converges to ip weakly in H 2 (SI). Then the sequence ipk converges 
strongly in Hq (St). We set y k := T (p> k ). Let v in JC (p) and set v k — sup (v, pf.) G K- (Vfe)- The sequence 
v k converges to v strongly in H^ (St). As y k := T (p k ), we get a (y k ,v k - y k ) > (/, v k - y k ), i.e. 

a (y k , yk) < a (y k , v k ) - (/, v k - y k ) . 

By Corollary [3]3] the sequence y k is bounded and weakly converges in Hq (SI) to some y. Using the lower 
semi-continuity of a, the previous relation gives 

a(y,y) < a(y,v) - (f,v-y). 
As y k > ipk, this implies that y > (p, therefore y := T (ip), □ 
We obtain the main result of this section. 

Theorem 3.4. Problem (V) admits at least one optimal solution p G H 2 (Sl) HHq(SI). 

Proof. Let <p k a minimizing sequence. As J (p k ) is bounded, p k is i? 2 -bounded and converges to <p weakly 
in H 2 (SI). By Theorem 13.31 the sequence y k := T (ipk) converges to y* ■= T (ip*) weakly in Hq (SI) and 
using the norms semi-continuity we obtain 

J(f*)-=h I (T &*) - z) 2 dx + % ( [ (Ap*) 2 dx) < liminfJ(^ fe ) = inf (V) . 
Jn \Jn J 

□ 

3.2 An Approximated problem (Vs) 

We use a trick of Barbu 0, and add adapted penalization terms to the approximated functional Js (here 
we add ^\\ip — ip*\\o) to force the relaxed obstacle family ip to converge to a desired solution p* of (V) . 
So for any S > 0, we define 



J s (ip) := \ 



(T s (p) -zfdx + v (J (Aipfdxj +\\ip-ip*\\ 



2 

L 2 (n) 



The approximated optimal control problem (V ) stands 

min{J 5 (^), peU ad }. (V s ) 

Theorem 3.5. Problem (V s ) admits at least one solution p s . Moreover, when 5 go to 0, the family p s 
converges to ip* weakly in H 2 (St), and y s := T s (<p s ) converges to y* ■= T (ip*), strongly in Hq (Si). 
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Proof. The functional Js is coercive, and lower semi-continuous on U. Therefore, the problem (V 6 
admits at least one solution p 5 . We set y s := T s ((p s ), and note, that for any S > 0, 



Js (^) < Js (<P*) ■= \ 



(T 5 (p*) - zf dx + v (J (&p>*fdx\ 



(3.1) 



By Theorem 13.21 we know that T 5 Up*) converges to T Up*) strongly in Hq (SI), so that Js (<p*) converges 
to J (p*) as S — > 0. Consequently, there exist So > and a constant j* , such that 

V 5 < S , Js (ip 5 ) < f < +oo. 

Consequently p 5 is if 2 -bounded uniformly, for any 5 < So- We use the Theorem 13. 2\ we get ip s converge 
to <p weakly in H 2 (SI) and strongly in Hq (Si) and y 5 converge to y := T (p) strongly in Hq (ft). As U a d 
is weakly closed, we have <p in U a d- By (|3.ip and the lower semi-continuity of Js, we get 



J(0) + i \\$-tp*\\l< liminfj 5 (/) 

d— >0 



< limsupj^ Up ) 

< limJx (p*) < limjfoj*) 

< J(®. 

This yields that \\<p — f*\\ < 0, then <p = p* and lim Js (<p s ) — J (p*) ■ 

5^0 

In addition this proves that any clusters points of p s is equal to p* , so that the whole family converges. 

□ 

3.3 Optimality conditions for problem (V s ) 

We give first necessary optimality conditions for problem (V s ) . Les us recall the following result on the 
Gateaux-derivative of the operator T 5 [1] . 

Lemma 3.1. The mapping T 5 is Gateaux-derivative at any p inU a d- 

V£ e H\ (Si) , T& ^ + T 0-' rS (jg) ^ ^ m H i {n) ; when T ^ Qi 

T 

where v 5 is the solution of the following equation 

Av s + /3' s (y s -p)v 5 = p' s (y s - <p) £ in SI, v s = on dSl. 
Proof. See pj. □ 
We define the approximate adjoint state p s in Hq (SI) as the solution of the following adjoint equation 

A*p s + Ps (y s - p) p s = y s - z in Si, P s = on dSl, 
where A* is the adjoint operator of A. As p s is the solution of the problem (V s ), we get 



d 
Jt 

That is 



VpeU ad , -rJs{^ + t(p-p s )) >Q. 



Vp G U a d, / (x S (y 5 -z)+ vAp s A (p - p s )) dx + (ip 5 - p*) (p - p s ) dx > 0, 
Jn Jn 
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where x 5 £ #o (^) an d satisfies 

A X 5 + & (y 5 - <p s ) x s = (v s - v>") (<p - <P B ) in fi. 

from the definition of p 5 , we obtain 

X 5 A*p 6 dx + [ ft {y 5 - tp 5 ) p s X S dx + v [ Atp 5 A {tp - <p s ) dx + [ {tp 5 - tp*) [<p - tp 5 ) dx > 0, 
q Jn J n Jn 

where A* denotes the adjoint operator of de A. Then 

A x s p s dx + [ f3'g {y 5 - if 5 ) p s x 5 dx + v [ A<p s A [<p - <p 5 ) dx + [ (<p s - tp*) (ip - tp 5 ) dx > 0, 
n Jn Jn Jn 

we obtain 

P's (y S - <p 5 ) P S (<P -<p 5 ) dx + v I Alp 5 A {tp - tp 5 ) dx+ I {tp 5 - tp*) {<p - tp 5 ) dx > 0. 
si Jn Jn 

In the sequel, we set 

/:=$(/-/)/eL 2 (0). (3.1) 

Finally, we obtain 

Theorem 3.6. Iff 5 is an optimal solution of (V 5 ) and y 5 :— T 5 {tp 5 ) , there exists p 5 in H 2 (Q) CiHq (f2) 
and p 5 in L 2 (Q) such that the following system holds 

Ay 5 + S {y 5 - p 5 ) = f in fi, y 5 = on dQ, (3.2a) 

A*p 5 + / = y 5 -z infl, p 5 = on dQ, (3.2b) 

(n s + tp s -(p;(p-<p s )+u(Atff,A(tp-(fP))>0, VtpEU ad . (3.2c) 

In the case U a d '■= L 2 (O), we make this optimality system more precise. 
Let x m W and choose (p = tp 5 ± x\ by the equation (|3.2c[) . we obtain 

{p 5 + <p 5 -tp*, X ) + is{Ap 5 ,A x ) =0, V X £W. (3.3) 
Set h 5 = Atp 5 in L 2 (ft), so that for any x m @ (O) , the relation (|3.3p gives 

{/i 5 + tp 5 — tp* ,x) + v (h? ') ^x) = ( m the distribution sens), 

that is 

-vAh 5 = pi 5 + tp 5 -tp* e%>' (Q) . 

Using the same techniques as in [6], we deduce that h 5 dn = 0. Consequently, h 5 S U, and it is the unique 
solution of 

-vAh 5 = f / + tp 5 - tp* in L 2 (ft) , /i 5 = on 90. 
The last relation may be written as 

—v {A 2 tp s , it) = (/i 5 , it) — (i^ 5 — <y9*, u) in O, (jC* 5 = on dfl. 
Thanks to Green's formula, the previous relation reads 
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— V 



A 2 ipudx — (fi s ,u) — ((p s — ip*, u\ — 

in 



So Aip vanishes on the boundary dfl, and we conclude that <p s belongs to W := {u\u 6 H 2 (ft) n 
Hq(Q,) et Au|aQ = 0}. Finally we have: 



Corollary 3.4. Assume conditions of Theorem \3.6\ are fulfilled. andU a d = U, then the optimality system 
(S s ) reads 

a (y s ,v) + (0 s (y s -<p 5 ),v) =(/,«), WeW, (3.4a) 
a* (p«, «;) + (//, w) = (y s -z,w), Vwe U, (3.4b) 

i/ (AVi «)-(/>«) = i VueW. (3.4c) 

-ffere a* denotes the adjoint form of a (associated with the adjoint operator A*). 

4 First order necessary optimality conditions for (V) 

In this section, we have to estimate p s , and gives more convergence results, then we may pass to the limit 
in the system l|3.4p as 8 — > 0. 

Theorem 4.1. When 8 goes to 0, p s converges to p* weakly in Hq (ft) and (/ converges to fi* weakly 
satr in H' 1 (ft) n M (ft), and 

(^*,P*)>0. 

where M. (ft) is the set of all regular signed measures in ft. 
Proof. Using l|3~4b"l) . we obtain 

a* ( P S , P S ) + I P' 5 (y S - <p s ) { P S ) 2 dx = (y s - z,p 5 ) . (4.1) 

As /?'(■)> 0, and thanks to hypothesis (|2"T2|) and lf!Q|) . we get 

ll^lljyi(O) - ^ lb ~ z ll_H"i(o) ' 



which implies that p s converges to p* weakly in Hq (ft). Consequently A*p s is uniformly bounded in 
H- 1 (ft) and 

fi s = -A*p s + y s - z. (4.2) 
Let 7 £ £ ^ x (M) be a family of smooth approximations to the sign function and satisfy the following [17] : 

7e (r) > W £ R, 

and 

{1 if r > e 
if r = 
-1 if r < —e, 

Then we can multiply (|4.2[) by 7 e (p 5 ) and integrate it over ft. As a result, we get 
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Letting e — > 0, we have 



// 7e (p s ) dx < C. 



\v ILi(n) - C '- 



Hence p 5 is bounded in L 1 (O) and consequently it is also bounded in jM (fi), thus, // converge to fi* 
weakly star in H^ 1 (O) nM(fi), such that 

+ ^* = y* - z in 0, p* = on 90. (4.3) 

As < (3' < 1, by using [|4.1|1 we get 

a* (p 5 ,p S )<(y 5 -z,p S ). 

And by the lower semi-continuity of a* 

(A*p*,p*) = a* (p*,p*) < liminf {y 5 - zj) = (y* - z,p*) . 

o — >oo 

From (|4.3p . we obtain 

< (A*p*,p*) = (y* - z,p*) - (fi*,p*) < (y* - z,p*) , 

so that 



□ 



<M>*) >o. 

In the sequel, we set £ s ■= (3s (y s — <p s ), then we obtain the following results. 

Theorem 4.2. When S goes to 0, £ 5 converges to £* weakly in L 2 (f2), where £* is negative and the state 
equation (|3.2a|) gives 

Ay* + r = /. 
Proof. From (|3.4a[l . we obtain with v = (3s (y S — f 5 ) 

{A (y s - <p*) , (3 S {y s - <p*)) + (fi s {y s - tp s ) , (3 5 (y s - <p')) = (/ - V, A (v* ~ /)) ■ 
For the seek of simplicity, we set r s := y s — ip s ; with l]2.ip . this gives 



+ E ^ * ( r ") ^ = (/ - V, ft (r 5 )) ■ (4.4) 

With hypothesis (H), we get 



8' x (r s )dx> I mV[^| 0' x (r s )dx. 

(4.5) 



> 0. 
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From P~H) . ((475]) and (H), we obtain 



|ft (r 5 )\\l 2m < ||/- VL 2(n) ||ft (^)L 2(n) +E IMi~(n) (H^^IL-c^) II& (^)IL 2( n)^ 
< max | 1,^ ||o,|| £oo(n) | (||/ - V|L 2(n) + ||Vr 5 || L2(a) ) (r 5 )|| L2(n) . 



Finally, we get 

lift (y* - || L2(0) < « (||/ - VL (n) + 1| v (/ - || L2(n) ) , 

where a := max / l, II jl°°(£7) ^ i so that 

II& - /)L 3(n) ^ a (11/ - V|| za(n) + 11/ - /Hhx ( o)) ■ 

Since (p 5 and y" 5 are respectively bounded in H 2 (ft) n (^) and Hq (O), we deduce that £ 5 is bounded 
in L 2 (fi), by passing to the limit where (5 — > 0, we obtain that £ s converge to £* weakly in L 2 (fi). Passing 
to the limit in (|3.2a|l , gives 

Ay* + r = /• 

where £*, is negative and we get y* G H 2 (O) (~l -Hq (^)- '-' 
Corollary 4.1. As W »" # 2 (^) H #o (O), 2/* : = ^* (<P*) &eZon 5 s io ff 2 (fi) n (Q). 
Proof. As C and / belongs to L 2 (ft), then Ay* 6 L 2 (ft) and y* G H 2 (ft). □ 
Now, we give some Lemmas, the proof the below Theorem 14,31 

Lemma 4.1. When 5 goes to 0, (/i 5 , (y 5 — i P S ) + ^j ^ (l 1 * \ V* ~ i P*)> an d (/■** > V* — </?*) = 0- 
Proof. By the definition of (3 and /i s l|3.ip , we get 

(//,(/-/) + )=0, 

where w + := max{0, v}, by Theorem 13.21 (y" 5 — ip s ) + converges strongly to (y* — (/?*) in TJq (TO; 
then 

(//,(/-/)+) ^( Ai *,y*-^*), 

and 



</i*,y*-^*} = 0. 

Lemma 4.2. Mien 5 yoes to 0, -> {C,P*), and (C,P*) = 0. 



□ 
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Proof. As befor, we set r s := y s — / , so that £ s = (3g (r s ). From the definition of (3 and (3', we get 
respectively 



= (Mr*) 



{-|<r«<0} 



r 4 )V<&; 



and 



(A*V) = (rY) = 



p s r 5 dx-2 / (r 5 ) /da; 

{^<-I} i{-I<^<o} V ' 



from that, we get 



Then, we obtain 



y{r*<-|} 



1/2 



As ff 1 (SI) ^ L« (0) with 2 < 9 < 6, we have 



j{r s <-^} ( 



1/2 



(/) 2 dx 1 < 



7 (p*) 4 



V? 



(?-2)/2g 



dx 



<C||/|| Hl(n) (measK<-I}) ( - 2)/29 . 



where meas {A} is the measure of the set A. Then we write 



(2 \ 1 I 2 
J i(r* + |) rfxj (meas{/<-I}) 



5 , 1 n (q-2)/2g 



We have 



1/2 



(4.6) 
(4.7) 



As (r 



lL 2 (fi) 



is bounded, we deduce that 



and 



£ / r 11 (r 5 + |) 2 dx<C, 

y{r«<-i} 

meas{r 6 < -±} < CS 2 . 
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then, by passing to the limit when S goes to 0, we have 

Jim (meas jr 5 < -|}) =0. 
Since Hp^Hh 1 ^) ls bounded, from (|4,6p and (|4.7p . we obtain 

| (£ S ,P S ) - | (mV) | < C (meas < - \ , 

and 

lim|(eV)-|(M 5 ,^)| = 0, 
so when 5 goes to 0, (£, s ,p s ) converges to 0, and with Lemma |4~2} we get 

<r,p*> = o. 

Finally, from the previous convergence results, we obtain the main result of this section 



□ 



Theorem 4.3. Let ip* , be an optimal solution of problem (P). Then Aip* belong to Hq (fl) and there 
exists p* in H% (Q), £* < in L 2 (ft) and fi* in H 1 (fl) n M. (fi), smc/i i/iai i/ie following optimality (S) 
system holds 



Ay* 


+ C = .f in ft, y* = on <9ft, 


(l.a 


A*p* + 


/i* = y* — z inVl, p* = on <9£1, 


(2.a 


z/AV - 


- n* = in Q, A<yS* = (p* = on 9fi, 


(3.a 




- </?*} = 0, 


(4.a 




<r,P*> = o, 


(5.a 




a*(p*,p*)-(z-y*,p*) <0, 


(6.a 




(p*,M*> >0. 


(7.a 



Remark 4.1. in /-?£]/ . Ito et Kunish had obtained the following optimality condition system (<S) 

Ay*+C = f, C =max(0,r+y*-^), (Lb 

A*p + /i* = y* - z in H^ 1 (ft) , (2.b 

(— vAip* + vip* + fi* , x — <£*) > for all \ £ U a d, (3.b 

(y* ~ ( < 3 *) = a - e - in ^, (4.b 

= a.e. in fl, (5.b 

a* (p*,p*) - (z - y*,p*) < 0, (6.a 

(n*,p*<f>)>0 for all (j> £ W 1 '' (fl) , with <f)>0, and q> n. (7.b 
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Indeed, they studied the following optimal control problem (V) 

(Ay-f,<f>-y)>0 for all cf>£lC(ip), 

with 

f £ L 2 (ft) , ip £ Uad with ip < in dfl, 
such that the cost functional is given by 

J{V) :=\\ (T(<p)-zfdx+%( [ (\p\ 2 + \Vp\ 2 )dx) , 

We notices a likeness between the two systems (S) and (S) , exepted for the equations (3. a) and (S.b) (are 
respectively the differential of the objective function J and J); i.e. in JWj , the authors treated the optimal 
control problem (V), such that ip belong to U a d '■= {p £ X : p (x) > 0, on dft and — a(ip,v) + (/, v) < 
(X,v) for all v € V with v > 0} with H 1 -obstacle (where A £ L 2 (ft) satisfying X > a.e. on ft), and 
in our work, we had study the optimal control problem (P), where ip is in U a d '■— H 2 (ft) n (ft), with 
H 2 -obstacle. 

Conclusion 

In this work, we treated the theoretical aspect of the problem (V), we proved the existence of optimal 
solutions and constructed a necessary optimality conditions system. Additional optimal obstacle regular- 
ity has been also provided. Currently we study the numerical aspect of the problem (V), via a numerical 
strategy based on the direct resolution of the optimality system and using a fixed point algorithm. 
The author is grateful to Prof. M. Bergounioux for their instructive suggestions. 
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